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Abstract 

The phase structure of hadronic matter at high density relevant to the 
physics of compact stars and relativistic heavy-ion collisions is studied in a 
low-energy effective quark theory. The relevant phases that figure are (1) 
chiral condensation, (2) diquark color condensation (color superconductiv- 
ity) and (3) induced Lorentz-symmetry breaking ("ISB"). For a reasonable 
strength for the effective four-Fermi current current interaction implied by 
the low energy effective quark theory for systems with a Fermi surface we 
find that the "ISB" phase sets in together with chiral symmetry restoration 
(with the vanishing quark condensate) at a moderate density while color su- 
perconductivity associated with scalar diquark condensation is pushed up to 
an asymptotic density. Consequently color superconductivity seems rather 
unlikely in heavy-ion collisions although it may play a role in compact stars. 
Lack of confinement in the model makes the result of this analysis only qual- 
itative but the hierarchy of the transitions we find seems to be quite robust. 
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1 Introduction 



The dynamics of nucleons in dilute environment at very low energy can be ap- 
proached via effective field theories such as chiral perturbation theory. For instance, 
the properties of two nucleons, both bound and scattering, can be very accurately 
calculated within the framework of well-defined strategy 0. The situation with 
denser systems such as heavy nuclei and nuclear matter is an entirely different is- 
sue. While the traditional phenomenological methods have been quite successful in 
correlating a wide range of phenomena, there is very little predictive power in these 
techniques in the density regime where new physics can set in. The question one 
would like to raise is what happens to a hadronic (nuclear) matter system when it 
is compressed to a high density as in relativistic heavy ion collisions and in neutron 
stars. This requires going, perhaps much, beyond normal nuclear matter density, 
Po ~ 0.16 fm -3 . At some high density, it is plausible that nuclear matter undergoes 
one or more phase transitions. For instance, neutron matter found in neutron stars 
could be unstable at a density a few times po against the condensation of kaons with 
an intriguing consequence on the structure of neutron stars and the formation of 
black holes [0, |3|]. At about the same density or higher, chiral symmetry restoration 
with a possible induced Lorentz symmetry breaking could take place 0, ||. 

Even more recently the old idea of color superconductivity proposed by Bailin and 
Love |J drew a revived attention in a modern context and became a focus of wide 
speculations and debates |7|, §]. If this phenomenon were to take place, it could 
do so in the regime where either kaon condensation and/or chiral phase transition 
and/or induced Lorentz symmetry breaking might set in and the question is how and 
in what language, macroscopic (effective) or microscopic (quark-gluon), adapted to 
QCD one should address this problem. So far the most economical way of describing 
dense hadronic matter up to and near the normal nuclear matter density is to treat 
the hadrons that enter as quasi-particles, that is, fully dressed particles, with certain 
symmetries (such as chiral symmetry and scale invariance) assumed to be relevant. 
This is the basis of Brown-Rho scaling || [TIJ which is found to provide one of 
the simplest explanations for the low-mass lepton pairs observed in the CERES 
experiments The CERES experiments probe densities a few times the normal 
matter density, so at least for these data, the hadronic quasi-particle picture seems 
to work. However in dealing with the phase transitions that go over to QCD's 
microscopic degrees of freedom, the hadronic quasi-particle picture should cede to 
one in which quasi-particles are quarks and gluons instead. The proper strategy 
must then be a Fermi-liquid theory of quarks and the phase transition of the color 
superconductivity type could be analyzed in terms of Landau Fermi-liquid theory 
and instability from a quark Fermi liquid. 
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Such an analysis involves the decimation of energy scales relative to the Fermi mo- 
mentum and a thorough study of the renormalization group flow of the couplings 
corresponding to operators which are marginal at the given energy scale [|12|]. If 
these coupling constants grow to order unity, one would have to resort to intu- 
ition to set up the operators corresponding to the new emerging physics |12|]. In 
the deconfmed phase of QCD at high density, the new physics which is controlled 
by one-gluon exchange is most likely the formation of a color superconducting gap 
reminiscent of diquark condensation. For a quantitative investigation of this effect 
see [0. On the other hand, we believe that the (possible) phase transition due to 



what we called "induced Lorentz symmetry breaking (or ISB in short)" in refs.[||, H 
occurs at rather low densities at which QCD is still strongly coupled. In this case, 
an educated guess of the relevant operators as required by the renormalization group 
analysis is difficult, and, if the induced symmetry breaking phenomenon interferes 
with quark confinement, impossible at the present stage. 

In this paper we shall therefore adopt a more coarse-grained point of view and ana- 
lyze the quark state at finite density resorting to a low-energy effective quark model. 
The model we shall consider consists of a gauge-invariant four-quark contact inter- 
action with the quantum numbers of one-gluon exchange supplemented in certain 
channels (e.g., the flavor-singlet vector channel) with collective modes associated 
with the presence of a Fermi sea. Such an interaction represents the lowest-dimension 
operator that would result when the gluon is integrated out from the QCD parti- 
tion function and would be dominant under the renormalization group flow as high 
energy states, both elementary and collective, are suitably "decimated." The con- 
straint we shall take into account in the current-current interaction - which controls 
the overall strength - is that it be consistent with the properties of pseudo-Goldstone 
bosons (pions) in matter-free space. The presence of Lorentz non-invariant collec- 
tive excitations could affect significantly the flavor- singlet vector channel (that we 
shall frequently refer to as Vo channel for reasons to be specified below) relevant to 
the ISB phase discussed in refs.||, |j. In particular, we shall take into account the 
possibility that the attraction in this channel could be enhanced by, and intimately 
correlated with, such collective modes as N*N~ 1 excitation considered in |L4[] for 
"dropping" vector masses in dense medium. 

We admit that the model is plagued with certain shortcomings (e.g. lack of quark 
confinement). Nonetheless we believe it yields a comprehensive view of all three 
possible phases, namely, chiral symmetry phase, the ISB phase and color supercon- 
ductivity phase. Given that lattice QCD measurements in the presence of chemical 
potential are totally lacking, qualitative studies of the sort we are presenting here 
seem to be the best one could do at the moment. 

The plan of this paper is as follows. We shall briefly review in section 2 the scenario of 



3 



chiral phase restoration at finite chemical potential with an emphasis of the induced 
breaking of Lorentz symmetry developed in The model we shall use is 

defined in section 3 wherein various condensates concerned are precisely defined. 
The effective gluon-exchange interaction that sets the order of magnitude of the 
parameters of the model is also given there. The phase structure that is expected to 
emerge as density and/or temperature is increased is described in section 4. Section 
5 contains discussions. Notations and definitions - some of which are quite standard 
- are given in the appendices to make our discussions self-contained. 



2 Phase transitions at finite chemical potential 



To render the paper self-contained and to clear up some of the confusion on the key 
points of refs.[|], 0], we briefly contrast the ISB scenario of chiral restoration^ to the 
more standard mechanism as provided by a constituent quark model. 



In the last decade, numerous investigations |L6| of constituent quark models have 
shown that the constituent quark mass M(p) rapidly decreases if the chemical po- 
tential n exceeds a certain value which we denote by fi c . In the constituent particle 
picture the density of states at a given momentum p, i.e. 

n(p) = 1 , E(p) = JMW + v 2 , (1) 

exp [(E(p) - n)/T\ + 1 v 

is governed by a sharp Fermi surface at zero temperature. In this case, the density 
of constituent quarks remains zero for \i < fi onset = M(/i). For M(/i) < ji < /i c , the 
density is a monotonic increasing function of /i. At the chiral phase transition at 
/i = fie, the constituent mass drops. An inspection of ([l]) shows that one expects a 
rapid increase of density. 

1 There seems to be some confusion caused by the term "ISB" which to the best of our knowledge 
has not been used before (apart from us) in the literature. Since we do not know any other way of 
calling it, we shall continue using this term although it could very well be a misnomer. The state 
of matter involved here was defined in ref.Q and clearly illustrated in a toy condensed-matter 
model in ||. One obvious question to ask is whether or not by "ISB" we are talking about a 
phase transition which is completely different from the chiral phase transition that is believed to 
be present in QCD. We should stress that the answer is no. That is, what we call "ISB" would 
be implicitly present if QCD were solved reliably (see e.g., fig. 2 in fll|). What we are doing here 
- which we believe is novel - is exploiting a characteristic of chiral phase transition that has not 
been identified before in the modeling of the QCD phase structure. As defined below, the "ISB" 
is characterized by the "jump" in baryon density as one goes across the critical chemical potential 
[i c corresponding to chiral phase transition and our point is that this is driven by an attraction in 
the vector channel caused by certain collective modes, in analogy to the potential introduced in 
the toy model of || (see figures 1 and 2). We will argue that this is related to the mechanism that 
may be responsible for "dropping" vector meson masses. 
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At finite density, our assertion is that collective modes due to vibrations of the Fermi 
surface - which must be present on a general ground when the renormalization 
group flow to the Fermi surface is properly studied in the effective Lagrangian or 



in hadronic models of collective excitations) 13 1 - must be taken into account for 
a correct description of finite density hadron matter. In 0,0, these modes were 
considered in terms of a local (effective) vector field V^, the vacuum expectation 
value (VEV) of the zeroth component of which measures the departure from the 
constant chemical potential /i and the discontinuity in the Fermion number density 
If the chemical potential reaches a certain strength fi ci , which is assumed to be 
smaller than /x c , an expectation value (Vo) forms and signals an ISB (on top of the 
explicit breaking via the chemical potential). In the constituent quark model, the 
density of states is given in this case by 

n(P) = exp[(E(p)-(V )-^/T} + l ■ (2) 

With the assumption (Vo) + > A*c the latter equation implies that the rise of 
the condensate (Vo) catalyzes the onset of the chiral phase transition. The state of 
matter after chiral restoration is, in the present case, considerably richer than what 
is given by the naive (standard) scenario: as explained in due to the dynamical 
origin of the condensate (Vo), light (non-relativistic) vector excitations are predicted 
to exist as pseudo-Goldstone-type modes which dominantly decay into dilepton. An 
intriguing possibility is that this light mode is just the lower branch of the two-level 
model of [[14[ . 



Due to asymptotic freedom, there are two important deviations in our effective 
quark theory from the naive constituent quark picture which was used above. Note 
that the low energy effective quark model which we shall employ below implements 
asymptotic freedom by a sharp 0(4)-invariant momentum cutoff. This analogy is 
of course only of qualitative nature and may not be quite realistic on a quantitative 
level. 

Firstly, the finite range of momentum integration implies a smeared Fermi surface 
even at zero temperature. Including a self-consistent treatment of the field Vq, small 
values of density might be observed also for /i < fi onse t- Technically, this rise of 
density as function of \i < {i onse t is a consequence of the smeared Fermi surface 0. 
From a physical point of view, this rise emerges from a polarization of the Dirac 



2 The behavior of density vs. chemical potential for a low chemical potential fi « Hd is 
irrelevant to the issue in question. For a dilute system, the correct degrees of freedom are nucleons 
or rather quasi-nucleons and not quarks or quasi-quarks. There is no reason to believe that an 
effective quark field theory without proper confinement mechanism can say anything about that 
regime which is of no concern to this work. 



5 



sea. Similar effects were observed before in NJL-type quark models when describing 
nucleons as solitonic excitations |T7 . 



Secondly, the next question one might raise is whether this quark system suffers 
from an instability signaled by a steady increase of (Vo). This instability is avoided 
by the shutdown of the quark interaction at large momentum transfer as expected 
from asymptotic freedom. This can be easily anticipated from the effective one-loop 
potential for Vq when one resorts to a quark model with sharp momentum cutoff A. 
At least for values Vq > A, one observes a monotonic increase of the potential as 
function of Vo telling us that the formation of large values of Vq is unlikely to take 
place. 

Finally, the quark interaction mediated by the "collective" field V^ is attractive (as 
suggested in |L4]) and thus favors the formation of the V condensate at sufficiently 



strong values of the chemical potential. One might question whether such an inter- 
action is not an ad hoc input foreign to the expected properties of low-energy QCD 
merely to induce an ISB mechanism. The answer to this question is this: firstly, 
at least for |Vo| < fi, the collective field describes perturbations of the Fermi sur- 
face. These fluctuations account for new important physics which comes into play 
at finite densities, and which is present besides the microscopic interaction mediated 



by gluons. At the hadronic level, the schematic model of |L4| suggests that these 
collective modes could be excitations of the A^*-hole type and other varieties - anal- 
ogous to phonons in solids and giant resonances in nuclei - in the isoscalar vector 
channel. The exchange of such collective modes between quarks could contribute an 
additional attraction. It may be possible to determine the size of the contribution 
from the collective modes using a model of the type studied in ref. [13]. But this 
program has not yet been fully worked out. In what follows, since we do not really 
know how to pin down the strength of the net attraction, we shall vary it within the 
range indicated by the one-gluon exchange plus that of collective modes and explore 
the qualitative structure of the possible phases. 



3 From the model to the quark ground state 
3.1 The model 

We shall use a path-integral technique, so we define our model in Euclidean space in 
order to ensure a proper convergence of the functional integrals. In particular, the 
quark fields will be Euclidean spinors. Our notation and conventions can be found 
in appendix |A|. The partition function which we will use throughout this paper is 
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given by 



Z 
L 

where 



J Vq Vq Vq c Vq c Va Vn a VA DA f £% exp I J d 4 x iA 
q \ fSo 1 kA^\ f q\ _ At A 



kA So 1 J 



9A- 



(3) 
(4) 



V u (x){-d' 2 8 tll/ + d ll d v )v v {x) - mlV^ 



9x 



(a - mf + 



7T 



q-1 



ip - a(x) + ij 5 7f(x)r + iV^(x)^ + ifMj 
Sq 1 := ip - a{x) + 27 5 7f(x)r T - iV^x)^ — ifi'y 



1 color j 
1 color 



(5) 
(6) 



where m is the current quark mass, \x the chemical potential and g&, g x and m v 
are (necessarily positive) coupling constants which control the strength of the four- 
quark interactions that result when the bosonic fields - that we shall simply refer 
to as "mesonic" fields - a, n, and A, At are integrated out. The role these 
couplings play will be discussed below. The parameter k e {l,i} will be chosen to 
reproduce the correct sign of the four-quark current-current interaction which is well 
known for a correct description of the low energy physics of the light mesons |T5| (see 
subsection |3.3Q. The charge conjugated quark fields possess the same kinetic energy 
as the quark fields and couple to the flavor singlet vector fields with opposite sign 0. 
The quark fields q, q and q c , q c are treated as independent fields. At the classical 
level, the equations of motion yield the familiar relations q = q* and q c = C(q^) T 
for Euclidean quark spinors (see appendix 0), where C is the charge conjugation 
matrix. 

Equation (|j) is the simplest Lagrangian that contains all relevant degrees of freedom 
for the physics we are interested in. The quarks form flavor doublets and belong to 
the fundamental representation of the SU(3) color group. The "mesonic" fields a, 
■K and Vfj, are color singlets, a and are flavor singlets while the pion fields form 
a flavor triplet. The quantum numbers of the A field will be given below. The 
fields a and i\ a are related to scalar and pionic degrees of freedom, respectively. 
The (collective) vector field describes vibrations of the Fermi surfaces and 
can be interpreted as the relativistic analog of the zero sound or the giant dipole 
resonance |J. Finally, A is the conjugate variable coupled to the diquark composite 
field and can therefore be viewed as the color superconductivity gap. 



3 As defined in (g) and (g), the collective field V M ~ (Vb,0) carries constants that reflect the 
possibly density-dependent) strength of the coupling. They are buried in the constant in 
One could think of the factor that we shall vary later reflecting, say, the 7V*A r ~ 1 -vector 
coupling strength of H. 
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Let us study the symmetry of the model @ in the case A, A^ = 0. In the chiral 
limit m — > 0, the Lagrangian L for the two light-quark system, (|j), is invariant under 
an SU(2)a x 5*?7(2)y flavor transformation of the quark fields q and q c , respectively 
(see appendix |C]). For a vanishing chemical potential fi, the model is 0(4) invariant 
corresponding to the Lorentz symmetry in Minkowski space. For /i ^ 0, the 0(4) 
symmetry is explicitly broken to a residual 0(3) invariance. 

In order to compare the meson-induced interaction of the model (HD with the low 
energy effective quark interactions, it is instructive to integrate out all meson fields. 
Defining 

A =: T5 (7) 

where 5 is a c-number and T specifies the Lorentz-, flavor- and color-structure of the 
superconducting gap, the following quark interactions (in derivative expansion and 
at low momentum transfer) are obtained 

+ mqq + ^-[qqqq - q^r a q q~f 5 T a q } - -^qYqqYq (8) 

+ mq c q c + ^—\q c qcMc ~ qclz>T a q c q c ^ 5 T a q c } - — — q c ^q c q c J^q c 
Ag x Am 2 v 

+ K 2 — q c TqqT^q c . (9) 
9A 

3.2 From the effective potential to the quark condensates 

The effective potential is a convenient tool for studying the ground state properties 
of an effective quark theory. Its global minima determine the condensates or ground- 
state expectation values of the bilinear quark fields that dictate the symmetries at 
low energies. The procedure is quite general and, in the case we are interested in, 
is equivalent to the BCS or Nambu-Gorkov formalism. 

Let us be more specific. Firstly we derive a non-linear meson theory by integrating 
out the quark fields. Note that there is an intrinsic relation between quark con- 
densates and ground state expectation values of meson fields. Exploiting the fact 
that the partition function @ is invariant under a shift of the mesonic integration 
variables yields the desired relations 

2g x (a-m) = (qq) + (q c q c ) 

2ml(V ) = (iqj°q) - (iq c J°q c ) 
-2g A (5) = (qcK-Tq) + (qnT ] q c ) . 

As expected, the quark and conjugate quark fields contribute to the scalar conden- 
sate with equal signs while they come with opposite signs in the case of the density. 



S 



(a) serves as an order parameter for the spontaneous breaking of chiral symmetry, 
while (5) can be directly interpreted as the gap of a color superconducting state. 
(Vo) is proportional to the quark density which is always non-zero for non- vanishing 
chemical potential. 

To define a convenient indicator - say, a "litmus paper" - visb which signals the 
onset of the ISB phase transition, we first note that, in the limit ml — > oo, the four- 
quark interaction in the vector channel vanishes, and is constrained to zero. The 
collective modes corresponding to the fluctuations of the field decouple. In this 
case, the density is entirely generated by the chemical potential. For finite values of 
ml, the interactions cause Vo to depart from zero. We define the indicator or signal 
for the ISB phase by 



PlSB 



p(V>,(,5>,<Vo>) - p(V),(<5),o) , (10) 



i.e. we subtract from the full ground state density p ((c), (6), (Vo)) the density 
p((a), (5),0) which is induced solely by the chemical potential, pisb therefore di- 
rectly measures the contribution of the collective dynamics to the density and is 
equal to zero for ml — > oo. While p ((cr), (S), 0) is expected to behave smoothly 
under increase of the chemical potential, Pisb summarizes the non-trivial behavior 
due to the ISB phase transition. 

To access the meson expectation values (in some approximation), one first introduces 
external sources that linearly couple to the meson fields in which we are interested. 
The generating functional in the case of our model (H) is therefore 

Z[j] = JVqVq Vq c Vq c Va Vn a VA Vtf £% exp j J d 4 x [L + L source ] 

Lsource = j a (x) a(x) +j (x)V (x) + j A (x) 5(x) . (12) 

The meson expectation value (Vo), and therefore the condensate (97 ?) will be non- 
trivial for non-vanishing values of the chemical potential p. For convenience, we 
introduce the compact notation j = (j a ,jo,jA) and = (a,Vo,S). The effective 
action is defined by a Legendre transformation with respect to the external fields 



A[4> c ] = - \nZ\j\ + J d 4 xj(x)0 c (x) , (13) 

where 

«*> - Wf ■ (i4) 

The effective potential U e ff(cj) c i) is obtained from ^4[0d] by resorting to space-time 
constant classical fields. 
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In the case of many quark degrees of freedom (e.g. large number of colors), the 
functional integral flTT|) can be approximately calculated neglecting mesonic fluc- 
tuations of the mesons (including A) around their mean-field values. Within this 
approximation, we find 

+ g x {a d -m) 2 + m 2 v v 2 d + g A 5 d . (15) 

It can be proven that the fermion determinant is real if the diquark vertex F is 
either hermitian or anti-hermitian. This will be always the case for the examples 
studied below. The classical fields at which the effective potential U e ff takes its 
global minimum yield the desired condensates. In the following, we will interpret 
a d m as quark condensate and v d m as quark density, and we will drop the index cl, 
since U e ff consists only of classical fields. We will study the phase structure of our 
model (|3|) by referring to different types of Lorentz- and color structures T of the 
diquark condensate. 



3.3 Gluon-induced quark interactions 

In this subsection, we will estimate the size of the parameters that figure in our model 
by comparing (§,[|) with the low-energy quark theory which may be obtained from 
effective one-gluon exchange. Note however that the resulting four-Fermi interaction 
represents a lot more physics than the exchange of one bare or dressed gluon and 
the effective size of the parameters will be fixed only when the additional corrections 
that arise in the decimation toward the Fermi surface (coming from, e.g., significant 
collective modes including the iV*iV" _1 excitation of |14j]) are suitably accounted for. 

Quark interactions which are mediated by the exchange of a dressed gluon can be 
described in terms of the action 



Sigiuon = - J d A x q [ip - Al T a J"]q- I d A x q c [ip + (T a f ^ 



Qc 



d 4 x d'y A#x) D£(x - y) A b u {y) , (16) 

where T a are the eight generators of the SU(3) color group. If V denotes the quark 
gluon vertex, the charge conjugate quark fields couple to the gluons via CT T C _1 0. 
Much work has been done on investigating the properties of light hadrons resorting 
to a wide variety of ansatze for the dressed effective gluon propagator D ab v [x — 
y) 0, |22], |3|]. Note that is assumed to be a positive definite operator to 
ensure convergence of the gluonic functional integral. When employed in the Dyson- 
Schwinger equation for ground state properties and in the Bethe-Salpeter equation 
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for description of light mesons, these effective interactions are known to provide a 
successful description of the spontaneous breakdown of chiral symmetry and the 
corresponding low-energy pion physics. Over the last decade, research has extended 
these effective 1-gluon-exchange models by using non-trivial quark-gluon vertices 
that follow from Taylor-Slavnov identities. These investigations have met with a 
remarkable success in low-energy hadron phenomenology |24|, |26|. 



Here, we shall not resort to such Dyson- Schwinger-type approaches to color super- 
conductivity, but will employ the corresponding low energy effective quark theory 
which emerges from an elimination of the high energy modes by means of a renor- 
malization group approach. In such an approach, the effective four-quark contact 
term provided by the current- current quark interaction naturally emerges, since this 
interaction respects all symmetries of low energy QCD. Moreover, it was observed in 
strongly coupled QED that such an interaction plays an important role in low energy 
physics due to a large anomalous dimension of the corresponding operator. In 
the context of low energy QCD, the current-current interaction dictates the Lorentz, 
color and flavor structure while its strength must be chosen to correctly describe the 
physics of the light mesons flnj. Such an effective low-energy quark theory coming 
from a renormalizable QCD-inspired quark model generically yields a correct phase 
structure although it may fail in the description of certain dynamical processes (for 
details see W> 



When applying the renormalization group decimation of energy scales to the model 
(|TH), the current- current interaction which is compatible with the low energy QCD 
symmetries is given by 

L = ...+ G[qlFT A rqql F T A ^q + q c lF(T A ) T ^q c q c l F {T A ) T Yq c 

- 2q c l F (T A ) T rqcqlFT A Yq}i (17) 

where its absolute strength G must be chosen for e.g. reproducing the correct pion 
physics in free space. 

We follow the ideas first presented in for the diquark phase, and assume that 
a diquark condensate respects chiral symmetry while it breaks the residual global 
color symmetry from SU(3) down to SU(2). The color, flavor and Dirac structure 
of the diquark vertices which we will investigate below are (see appendix |B|) 



r 




(scalar) 


r 


= e lk e^ 


(pseudo-scalar) 


r 


= e ifc e a/33 7 M 75 


(vector) 


r 


= e ifc e a/3 V 


(axial- vector), 



where i, k — 1, 2 act in flavor space, and a, (3 = 1 ... 3 are color indices. Equation 
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(|17|) will be our starting point for estimating the coupling strengths of our model 
(^]). For this purpose, we rearrange the four-quark contact interaction ( pi) by means 
of a Fierz transformation so as to make sure that the classical equations of motion, 
i.e. the mean-field level, of our model equal the Dyson- Schwinger equations with the 



interaction (|17|) . We relegate the details of this calculation to appendix p.l| . The 
result is 

G 

L = ...+ — (ql F l c q ql F l c q - q l c -f 5 r a q q l c75 r a q ) (18) 



G 



1 1 

2N ~ 8 



G 



qlplc^q qlplc'fq 
{q )7 e ik e^^q P k qfe lm e^ ( & ) 



AN 

+ §j (tit e ik e^ q? q? e lm e^ {q c f m 

- §j (&)? e»c^ 75 gf gf e im e^ 75 , 

where 1 F , 1 c are unit matrices in flavor and color spaces, respectively, and G sets 
the scale of the interaction strength. N = 3 is the number of colors. Comparing 
this result with ([5],^), we first deduce the "sign" of the meson quark couplings, i.e. 

^scalar ^ > ^ps 1 > ^vector i j ^axial 1 i (^^) 

and then read off the corresponding strengths, i.e. 

1 2 N 1 

9 * = 2G ' m ^ = Gi^vyi ' (20) 

scalar ps _ ~ vector axial _ 4iV 

9a — 9a — ~q > Pa — 9a — ~q~ ■ l z -U 

In particular, we find the following ratios for N = 3 

g x : m 2 v : ^a^ = Pa^ = 1:6:6:12. (22) 

In the next section, we will study either case of diquark condensation and will tacitly 
assume that a mixed phase, e.g. with the scalar and vector diquark condensates 
simultaneously present, will not occur. 



The parameter ratios fl22|) represent the quark interactions carrying the quantum 
numbers of one gluon exchange. As explained in section |2|, the attraction already 
present in the current-current interaction for the Vq channel with three colors and 
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two flavors (see Appendix D for a caveat that may follow from a different Fierz rear- 
rangement) is likely to be further enhanced by collective modes. In our calculation, 
this effect will be duly accounted for by lowering the mass m^. Since little is known 
at present about the size of this additional attraction, we will explore in the next 
section the phase structure for a wide span of parameters while the order of mag- 
nitude of the parameters will be kept at the values given by the effective one-gluon 
exchange. 

4 Phase structure 

At finite baryon density, a rich phase structure is expected to appear in QCD. 
Whereas at zero density a quark condensate is formed in the QCD quark sector, two 
possible phases might exist at finite values of the chemical potential: the ISB phase 
in which Lorentz symmetry is "spontaneously" broken in the sense defined in ||, H 
and the phase in which a diquark condensate is formed which breaks the residual 
global color symmetry of QCD 0. In the following subsections, we will explore the 
phase structure which arises from our low-energy effective quark theory when the 
chemical potential is increased. 

4.1 Scalar and pseudo-scalar diquark condensations 

We follow the ideas first presented in J7| for the diquark phase. There, it was 
conjectured that the color superconductor instability occurs in the scalar diquark. 
In this section, we will extend the analysis of to the study of possible scalar, 
pseudo-scalar, vector and axial-vector diquark condensations as well. 

In this subsection, we will study the possible emergence of scalar as well as pseudo- 
scalar diquark condensates, i.e. 



Here i,k = 1 ... 2 are flavor indices, while a, (3 = 1 ... 3 are color indices. Dirac 
indices are not explicitly shown. The behavior of the Euclidean quark spinors under 
Lorentz transformations can be found in appendix [B]. If one of these condensates is 
realized in a finite-density quark matter while the quark condensate is zero, chiral 
symmetry is restored (see appendix |C.2| ). 

We first consider the case of the pseudo-scalar superconducting gap, A = €ik€ al33 1 5 
in conjunction with the ISB phase transition extensively discussed in [|[ |j. The 




(scalar) 
(pseudo-scalar) 



(23) 
(24) 
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fermion determinant in (]i~5f) can be calculated in a closed form (see appendix |E|), 
i.e. 



Det ps = J] (k 2 + + a 2 ) U (k 2 _ + a 2 ^ 

it 

[(w + A; - z(/i + f ) w-) 2 + w+A; 2 + w 2 cr 2 § (25) 



where 



5 



2 



fc-t := (k ± + v) , k) T , w± := 1 ± 

The effective potential is 

Uff f (a,8,v) = - In Det ps 

+ g x {a-m) 2 + m 2 v v 2 + g A 5 2 . (26) 

The model (^) is considered to be an effective low-energy quark theory in that the 
momentum integration implicit in ( f26|) is cut off at the scale A. Note that for a 
reliable investigation of the model at finite values of the chemical potential it is 
crucial that the momentum cutoff leave the Lorentz, i.e. 0(4), symmetry intact. To 
assure this, we used a sharp 0(4) invariant momentum cutoff. 

In the case of the scalar diquark composite, we choose A = ^ike" 133 5. The func- 
tional determinant in ([T^) is given by (see appendix [E]) 

Det scaiar = J] ( k l + ° 2 ) 12 (27) 

k 

{k 2 _ + a 2 ) (w + k — i(fj, + v) wJ) 2 + w 2 + k 2 + w 2 + a 2 



Comparing Q27D with ( f23| ) , one immediately sees that the determinants with pseudo- 
scalar and scalar diquark entries, respectively, coincide if chiral symmetry is re- 
spected, i.e. cr = 0. In the following, we will assume that we will not encounter a 
strong CP-problem, but that the scalar diquark condensation takes place. 

We explored the parameter space, looking for the global minimum of the effective 
potential U e ff, by resorting to Monte Carlo techniques. We find generically three 
different minima of the effective potential. One minimum shows a strong chiral 
condensate while the superconducting gap and v are close to zero. We call this 
minimum min — a. The second minimum, so-called min — v, is located at strong 
values of v with the order parameters cr, 5 being small. Finally, the third minimum 
is dominated by a large value of the superconducting gap 5, and we shall call this 
min — A. All these minima are possible candidates for the true ground state which 
is characterized by the global minimum of U e ff. 
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Figure 1: Scalar diquark case: the effective potential U e ff in units of A 4 /47r 2 of 
the local minima as function of the chemical potential \x for m\ — > oo (left) and 
m l = 9 X = 9 a (right). 



To establish the instability due to color superconductivity at high density as pro- 
posed in 0, we turn off the interaction strength in the vector channel by tuning 
m 2 , — > oo and study the competition between the chiral phase and the superconduc- 
tor phase as a function of the chemical potential. The numerical result of our model 
is shown in figure p] (left panel). It indicates that the scalar diquark condensation 
sets in at some large values of the chemical potential fi = 0(A), where the use of 
a four-Fermi contact interaction is most likely unjustified. Since we are here only 
interested in the qualitative phase structure, we shall not consider the non-local 
interaction due to the full effective one-gluon exchange. 

The situation becomes more dramatic if we use a moderate value for m\. For instance 
for gA = 9x = m 1 = A 2 /47r 2 , the effective potential at the minima comes out to be 
as shown in figure |I| (right panel). It turns out that ISB phase transition takes place 
at a small value (/i ~ 0.4A) of the chemical potential, whereas the minimum which 
is characterized by the scalar diquark condensate, i.e. rain — A, turns into a saddle 
point. 
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Figure 2: The chiral condensate a, ISB indicator pisb and superconducting gap S as 
function of the chemical potential for the 1-gluon exchange (left) and an unnatural 
weak four-quark interaction in the vector channel (right). 



In the case of finite m 2 , we resorted to a full Monte- Carlo treatment to calculate the 
strength of the order parameters a, v and 5 vec . It turns out that the scalar diquark 
condensate is not important for ground state properties if we choose the typical 
interaction strength of order unity, i.e. g x = g& = ml = A 2 /An 2 (see figure 0). 
At zero chemical potential, the quark condensate (a ^ 0) signals the spontaneous 
breakdown of chiral symmetry. The ISB phase (min — v) and the superconducting 
phase (min — A) appear as meta-stable states. When the chemical potential is 
increased above a certain critical value /x, a first-order phase transition takes place 
from the chiral broken phase to the ISB phase. An "induced" breaking of Lorentz 
symmetry takes place while chiral symmetry gets restored. This mechanism as well 
as its consequences for the light particle spectrum were extensively discussed in 0, § . 

We then reduced the strength of the quark interactions in the vector channel by a 
factor of 4, i.e. g x = <7a = A 2 /47r 2 , ml = 4A 2 /47r 2 , as this is the channel that gives 
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Figure 3: The effective potential U e ff in units of A a /4ti 2 of the local minima as 
function of the chemical potential /i for m\ — > oo. 



rise to the ISB phase. In fact, we searched for a region of parameter space where the 
superconducting state constitutes the state with the lowest energy density. In this 
case, one indeed observes a superconducting state at large values of the chemical 
potential (see figure 0). It appears, however, that if one takes the interaction strength 
to be of order unity - which is a typical strength one expects as before, the color 
superconductivity is unlikely to play a role in the density regime relevant for the 
physics we are interested in. This is because the ISB phase has a lower vacuum 
energy density. 

4.2 Vector diquark condensation 

One might wonder whether a vector diquark condensate and/or an axial- vector 
diquark condensate, i.e. 

(ggi^VTs^) , (vector) (28) 
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(qci ^ e a/33 l^ql ) , (axial-vector) , (29) 

might not compete with the ISB phase transition if the interaction is given by the 
effective one-gluon exchange. Since both diquark condensates transform as vectors, 
one might suspect that they could be more sensitive to the chemical potential which 
transforms as the fourth component of the vector. 

The fermion determinant in (|15|) can be readily calculated for the vector diquark 
(A = €ik€ af33 7 75 5) and the axial- vector (A = e^e"^ 3 7 75 5) diquark (see appendix 

@), 

DeW = I] { k l + ° 2 ) 12 (k 2 _ + cr 2 y (30) 

k 

Uw-ko — + v) w + ) 2 + w 2 + k 2 + w 2 _a 2 
Vet axial = Y[ (k 2 + + a 2 ) 12 (k 2 _ + a 2 ) A (31) 

k 

(w-k — + v) w + ) 2 + w 2 + k 2 + w 2 + a 2 

In order to gain some insight into the competition between the chiral phase and the 
phases exhibiting the various diquark condensates, we studied the limit m 2 — > 00 
which amounts to the suppression of the ISB mode. The effective potential of the 
various local minima - which are the possible candidates for the true ground state - 
is shown in figure |3| as function of the chemical potential. It turns out that the scalar 
diquark condensation sets in at [A — 0(A) whereas the other diquark condensates 
merely correspond to meta-stable states. 

The reason why the vector diquark condensates have no significant influence on 
the ground state properties is that the four-quark interaction in the vector diquark 
channel has the "wrong" sign. If for the purpose of illustration one changes this 
sign by setting K 2 ectOT = 1 (compare with (|19"|)), the vector diquark order parameter 
has a much stronger effect (see figure ^), but still cannot prevent diquarks from 
condensing in the scalar channel. 

4.3 Finite temperatures 

Finally, we study the impact of temperature on the magnitude of the order pa- 
rameters a, A and the ISB indicator pisb- In view of the results of the previous 
subsections, we will only focus on the emergence of the scalar diquark condensate. 

In taking into account temperature effects, we use the imaginary time formalism by 
introducing Matsubara frequencies and by replacing the ko momentum integration 
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Figure 4: Phase diagram for a medium size m 2 = A 2 /4n 2 interaction strength in 
the vector channel (left) and for a small size interaction m 2 = 8A 2 /47r 2 (right). 



in (|30|) by a discrete sum, i.e. 

k = (2n + l)7rT, f dk -> 2nT £ , (32) 

where T is the temperature. Note that the momentum integration is truncated 
by an 0(4) symmetric momentum cutoff A, i.e. k\ + k 2 < A 2 . This implies that 
the momentum integral vanishes for 7rT > A. From a physical point of view, this 
amounts to assuming that the non-local (e.g. gluon-induced ) interactions at large 
momentum transfer (k 2 > A 2 ) that match the four-Fermi interaction of our effective 
quark theory at the scale A contributes only perturbatively to the order parameters 
a, v and A (see []i"9"| ). We therefore expect the melting of all condensates at T c « 
A/7T. Using a generic energy scale of A m 600 MeV, one finds roughly T c m 200 MeV 
(at zero density), which is the right order of magnitude as revealed by lattice QCD 
studies [ pC|| . 

Using Monte Carlo techniques, we searched for a given chemical potential and a 
given temperature for the global minimum of the effective potential U e ff which 
occurs for the parameters a, v and 6. We then calculated pisb ( |10D for the signal of 
the ISB phase transition. Even at a medium scale of the chemical potential (// < A) 
one observes a significant jump of Pisb- Although, strictly speaking, there is no 
phase transition but a crossover, one can define rather sharp phase boundaries in 
parameter space by 

c > Vjsb, cr > 5 chiral phase 

visb > c, visb > ^ ISB phase 

5 > a, 5 > visb color superconductor phase, 

where visb '■= Pisb/^^1- The numerical result of our model is sketched in figure |j. 
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One finds that for the interaction strengths g x , g^ alar ; rn% all of order unity, diquark 
condensation plays no role. It is the ISB phase Pisb 7^ that characterizes the quark 
state. At large values of the chemical potential, we expect the contribution from 
the quark interactions to the density to become small, and that there is a smooth 
transition from the ISB phase at some medium values of the chemical potential to 
the phase of perturbative QCD where the density is largely induced by the chemical 
potential. 

If one were to choose an unnaturally weak vector interaction strength 
m l ^ 9x^9A alar i one m ight observe "shells" built up from the ISB phase and the 
color superconducting phase, respectively. We note, however, that the parameter 
range \x ~ A is beyond the scope of the present model. In order to access these ranges 
of the chemical potential, a non-local and asymptotically free quark interaction (as 
proposed e.g. in [OH]) might be required. 



5 Discussions 

A detailed analysis of the possible phase structure of finite density quark matter 
based on an effective low-energy quark theory is presented. The four-quark contact 
interaction used in the study is provided by the effective quark current-current in- 
teraction which naturally emerges from a renormalization group decimation of high 
energy scales in a many-body system with a Fermi sea. From a phenomenological 
point of view, the effective current- current interaction is well known for a correct 
description of the physics of the light mesons at least at low densities. The simplest 
way to understand what we have done in this paper is that our strategy corresponds 
to interpolating the effective bottom-up partition function (|) with the explicit col- 
lective fields and A (as well as the pseudoscalars tt and the scalar a) to that 
with the four-Fermi interactions (|T7D that must arise from a top-down approach 
starting from the QCD partition function. Although the effective theory might be 
incomplete in some aspects of nonperturbative structure of QCD, such as lack of 
quark confinement and additional many-body effects, we believe that it qualitatively 
describes the phases of quark matter when the chemical potential is varied. Hope- 
fully such an approximate description of quark matter phases at high density will 
serve as an important input for a thorough renormalization group study (based on 
a decimation of energy scales relative to the Fermi momentum [TJJ) which will then 
yield more quantitative results. 

Our finding is that the generic interaction strength provided by the current- current 
interaction in the diquark channel is much too weak to trigger a scalar diquark 
condensate. Pseudo-scalar, vector and axial-vector diquark condensates are even 
less favored. Instead of the color superconductivity, we find an induced Lorentz 
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symmetry breaking to occur at medium values of the chemical potential leading to 
chiral symmetry restoration with color symmetry remaining unbroken. 



In sum, we are led to the following conjecture for the quark phase structure: at 
zero density, strong QCD interactions break chiral symmetry spontaneously leading 
to quark condensation. At small finite values of the chemical potential, the quark 
condensation weakens in the presence of Lorentz-symmetry-breaking terms, and a 
"roll-over" from (qq) to (q-foq) is observed (called "ISB phase transition"). One 
of our main points is that this phenomenon in quark language is related to the 
hadronic description of the chiral phase transition and the "dropping mass" of the u 
meson in JTJJ. As the chemical potential is further increased, the quark interaction 



is weakened due to asymptotic freedom. At a certain (second) critical value of 
the chemical potential, the interaction is no longer strong enough to sustain quark 
condensation in one of the Lorentz channels. In this case, it is most likely that color 
superconductivity occurs since an arbitrarily weak (but attractive) interaction at 
the Fermi surface is sufficient for (color) Cooper pair formation. 

We should mention that our description may be much too simplisitic to be realistic. 
It has been shown by Ilieva and Thirring |28| that even in a simple quantum me- 
chanical system involving both the ISB-like (called "mean field" by the authors of 
28f] ) phase and superconducting phase, the phase structure can be quite intricate 



and complicated. It is perhaps significant to note in particular that as in the Ilieva- 
Thirring model, an attractive "mean field" potential could even trigger (instead of 
prevent) superconductivity in contrast to what we have found in this paper. It is fair 
to say that the situation in QCD could be immensely more intricate and complex. 

Finally let us comment on the apparent breaking of the local gauge invariance by the 
formation of the diquark condensate. In a recent important work, Schaden et al. |2"9 ] 



exploited a topological field theory to obtain a complete gauge fixing while avoiding 
the Gribov problem. In the so-called equivariant gauge, gauge invariance was fixed 
up to a global symmetry in color space. The authors then showed that quark con- 
finement (in these gauges only) could be induced by a spontaneous breakdown of this 
residual global color symmetry. As a consequence, owing to Goldstone's theorem, 
light colored fields (which arise in these gauges) decouple from local gauge invariant 
operators but contribute to the Wilson loop (due to its non-local nature) |29| . These 
Goldstone bosons then could produce an area law (and therefore confinement) due 
to inherent infra-red singularities 



Except on a lattice, any approach to color superconductivity, be that via an effective 
quark theory or via an renormalization group analysis, will have to be based on a 
completely gauge fixed QCD. In this respect, diquark condensation does not break 
local gauge invariance, but breaks instead the residual global color symmetry. In 
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view of the results presented in [29|, one encounters the interesting possibility that 
a hypothetical phase transition to a color superconducting state at high density 
restores chiral symmetry, but does not imply quark deconfinement . However, we 
should recall that, at least in certain supersymmetric Yang-Mills theories, chiral 
restoration and deconfinement occur simultaneously 
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A Notation and conventions 

The metric tensor in Minkowski space is 

g„ v = diag(l, -1,-1,-1) . (A.l) 
We define Euclidean tensors T(E) from the tensors in Minkowski space T(m) by 

t$*" = « r hy Tfc*« Ul ___ un , (a.2) 

where r and s are the numbers of zeros within . . . fi^} and {v\ . . . u n }, respec- 
tively. In particular, we have for the Euclidean time and the Euclidean metric 

X °(E) = ix °(M), 9(E) = diag(-l, -1,-1,-1) . (A.3) 

Covariant and contra- variant vectors in Euclidean space differ by an overall sign. For 
a consistent treatment of the symmetries, one is forced to consider the 7^ matrices 
as vectors. Therefore, one is naturally led to anti-hermitian Euclidean matrices via 

(H), 

7(% = H( M ) , 1(E) = 7(M) • (A.4) 

In particular, one finds 

(7fe) f = -7&), {7fe,7(V = He) = "2<W- (A.5) 

The so-called Wick rotation is performed by considering the Euclidean tensors (|A.2[) 
as real fields. 
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In addition, we define the square of an Euclidean vector field, e.g. V^, by 

V 2 := = ■ (A.6) 

This implies that V 2 is always a positive quantity (after the wick rotation to Eu- 
clidean space). 

The Euclidean action Se is defined from the action Sm in Minkowski space by 

exp{iS M } = exp{S E } ■ (A.7) 



Using QA.2|) , it is obvious that the Euclidean Lagrangian Le is obtained from the 
Lagrangian Lm in Minkowski space by replacing the fields in Minkowski space by 
Euclidean fields, i.e. Le = Lm- 

Let the tensor A^ v denote a Lorentz transformation in Minkowski space, i.e. 

AVV> Q/3 = <T ■ (A.8) 

Using the definition (|A.2|) , one easily verifies that Mf E \ v are elements of an 0(4) 
group, i.e. Aj E ^A( E ) = 1) which is the counterpart of the Lorentz group in Euclidean 
space. 

In order to define the Euclidean quark fields, we exploit the spinor transformation 
of the Euclidean quark field 

9(B) (x E ) -»• q'(E) { x 'e) = S(A { E))q(E)(x E ) , x E ^x' E = Ax E , (A. 9) 
S(A iE) )^ E) S^A iE) ) = (A^Tfo, (A.10) 

where the matrices are given by 

i 

S(A {E )) = exp{--u flu a^ ) } (A.ll) 

for det(A) = 1 only. It is obvious that one must interpret 

Q(E) = q\ E) (A. 12) 

in order to ensure that e.g. the quantity <?(£) 7^ <?(£:) transform as an Euclidean 
vector. We suppress the index E throughout the paper and mark tensors with an 
index M, if they are Minkowskian. 
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B Charge conjugation in Euclidean space 

The equation of motion of a Dirac spinor which couples to a U(l) gauge field is 
given by 

-id^q\x)Y - q^(x)e4 = 0. (B.l) 
Defining the charge conjugated spinor by 

q c (x) := C (q\x)) T (B.2) 



with some Dirac matrix C which will be specified below, equation (|B . 1| ) can be cast 
into 

CifC~ l q c {x) + eC4 T C~ l q c {x) = 0. (B.3) 

One finds that q c (x) satisfies a Dirac equation with reversed charge e — > — e if C 
satisfies C7JC -1 = —7^- Let us choose a standard representation of the Dirac 
matrices, i.e. 

with 77 the SU(2) Pauli matrices. Note that 

{7^,75} = 0, 7 t = J? = 75 . (B.5) 

Using this representation, one finds 

C = 7° 7 2 , & = C- 1 = C T = -C . (B.6) 

Let us investigate the behavior of the charge conjugated spinor q c under Lorentz, 
i.e. 0(4), transformations ( |A.11[ ). Replacing S(A)q for q in the Dirac equation (B.l) 
and re-deriving q c from this equation, one finds that the charge conjugated spinors 
transform as 

q' c = C (S-\A)) T C-'qc. (B.7) 

We now confine ourselves to Lorentz transformations ( |A.11| ) with det(A) = 1 and 
will check the behavior of q c under parity transformations later. Firstly note that 

{l^lu - lvln) T = llll - llll = ~ C {lul^ ~ 7^) C' 1 , 
and therefore 

S T (A) = CS\A)C~\ (detA = l). (B.8) 

Using the last equation, one readily verifies for a solution q c of the equation of motion 
that 

q' c = C(gt<?t(A)) = c(&) = S(A)q c . 
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An inspection of (|A. 10| ) shows that 5(A) = 70 generates a reflection of spatial 
coordinates and therefore serves as an example for a Lorentz transformation with 
detA = — 1. We learn from ( p.7| ) that the conjugated spinor behaves as q c — ► — q c 
under this reflection transformation. We therefore conclude that the quark spinor 
q c transforms as 

q' c = det(A)S(A)g c . (B.9) 

In particular, one finds that the diquark condensates (<? c 75 ( z) an d {qd^l^Q) transform 
as a scalar and a vector, respectively, as in Minkowski-space. 



C Chiral symmetry in Euclidean space 



C.l Generalities 

For simplicity, we will illustrate the emergence of chiral symmetry in effective meson 
theories for the case of two quark flavors. Generalization to Nf quark flavors is 
obvious. One observes that the quark kinetic energy, i.e. 

q(ip - a + i75vrr) q (C.l) 

is invariant under a simultaneous transformation of meson and quark fields. Let V 
be an SU(2) matrix which acts only on flavor indices, and V a matrix which exists 
in flavor and Dirac space, i.e. 

V = exp {iar} , V = exp {iarlo} , (C.2) 

For convenience, we introduce the flavor matrix 

M = a — inr . (C.3) 

Then the invariance under V G SUy(2) transformation is observed for 

q' = Vq, q = qV ] , M' = VMV ] . (C.4) 

To display the axial- vector symmetry, we choose U to be another SU(2) matrix 
acting in flavor space and define 

U = exp{i(3r} , U = exp{i/3r 75 } . (C.5) 

Invariance under axial- vector chiral rotations, U G SUa(2), is generated by 

q > = Uq, q = qU , W = UMU . (C.6) 
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The so-called chiral circle is invariant separately under both transformations, i.e. 

a 2 + vf 2 = itr (MM f ) . (C.7) 

Note, that q and q are independent fields. This makes it possible to assign different 
flavor matrices to either field. Equations (|C.4j ) and (|C.6]) establish the well known 
SUa(2) x SUv(2) chiral symmetry. 

There is a subtle difference between the Euclidean formulation and the formulation 
in Minkowski space. Let ip and tp denote the spinor solutions of the equations of 
motion. In Minkowski space, one finds 

= ^(m)7(°m) , ${m) = (W(A/)) f 7(V) = $(M) u ( c - 8 ) 

Equation (|C.8| ) implies that the saddle points are degenerate by axial chiral rota- 
tions. This situation is different in Euclidean space. There one finds 

= ^\e) > $'{E) = ( w ^(s)) f ^ u ( C - 9 ) 

implying that the degeneracy of the saddle points is lifted. Since we are interested 
in the non-linear meson theory that results from integrating out the quark fields, 
the resulting fermion determinant depends only on the chiral circle ( |C.7| ) in either 
formulation. 

To illustrate this point for the Euclidean case, we calculate the fermionic determinant 
for the space-time constant fields a, n. The eigenvalues of the Dirac operator in 
momentum space, i.e. Ijt — o + mf 75, are given by A± = — o ± i\Jk 2 + tt 2 . One 
therefore obtains the desired result 

Bet = U(\ + \-) 2 = Y[(k 2 + a 2 + n 2 ) 2 . (C.IO) 

k k 

Finally, we must specify the behavior of the other independent quark fields q c , q c 
under chiral rotations. We define 

Qc = (V f )% c , t = QcV T VESU V (2), (C.ll) 
q' c = U T q CJ q' c = qM T U G SU A (2) (C.12) 

Using (|C.4|) and (|C.6|) , one readily shows that the kinetic energy for the charge 
conjugate quark fields, i.e. q c S^ 1 q c , is chiral invariant. 



26 



C.2 Chiral behavior of diquark fields 

In this subsection, we study the behavior of the composite field q c e T q under chiral 
transformations (|C4| , |C6| ) and ( |0 . 1 1| J07T2|) , respectively, where T is a Dirac matrix 
which will be specified when needed, and is the totally anti-symmetric tensor 
acting in flavor space only. Note that here q, q, q c , q c are treated as independent 
fields with the transformation properties as given above, and that one must not use 
the classical equation of motion, which e.g. yields the relation q c = C(q*) T . In 
the Minkowskian formulation there would be no difference, since chiral symmetry is 
manifest at classical level. As explained above, this is not the case for the Euclidean 
formulation. 

It is convenient to introduce left- and right-handed projectors by 

Pl = ^(1 + 75), P*=^(l- 7 6), PI = Pl, Pr=Pr- 

(C.13) 

The chiral matrices V G SU V (2) Q and U G SU A (2) O can be written as 

V = V x l D , U = U x P L + W x P R , (C.14) 
where U and V are matrices with iso-spin indices only. 
Let us study chiral vector transformations first. In this case, we find 

(q c erq)' = q c V T eTVq = q c eT q , (C.15) 

where we have used 

tikVilVkm = Qm • (C.16) 

The diquark composite is invariant under chiral vector transformations. 
In the case of chiral axial-vector transformations, one obtains 

(q c eTq)' = q c U T eTUq = q c (u T P L + (wf P fl ) eT (lJP L + U*P R ) q. 

. (C ' 17) 

In the case of an axial (T = 1) or scalar (T = 75) diquark composite, we find that 
the axial-vector chiral symmetry is respected, i.e. 

(q c eTq)' = q c U T eUP L q + q c (ll^f eWP R q. (C.18) 

Using ( |C.16| ), the invariance under axial- vector chiral transformations becomes obvi- 
ous. In the case of an axial- vector (T = 7^), the diquark composite is not invariant 
under chiral axial-vector transformations, i.e. 

(Qc q)' = q c U T e C/ V Pr Q + Qc (^f e 7 M UP L q . (C.19) 
The same is true for vector diquark composites (r = 7^75). 
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D Collective field approach to quark dynamics 



D.l Fierz transformation of the current-current interaction 



As emphasized in numerous places in the main text, the dynamics in the Vq chan- 
nel associated with the ISB effect is not dominantly controlled by the exchange of 
one dressed gluon. However in the literature, analysis of the QCD phase structure 
involving color superconductivity are made with one-gluon exchange four-Fermi in- 
teractions. In this appendix, we write down the various terms that result when the 
current-current interaction coming from one-gluon exchange is Fierzed. In doing 
this, we shall make use of the relations 



rpA rpA 

1 af3 1 j8 



$ik Sim 



1 

N 



N Oik Oim + N T lk T im , 



(D.l) 
(D.2) 



where T A are SU(N = 3) color matrices with the proper normalization, i.e. 
tr{T A T B } = 5 AB /2, and r a are the SU(Nf) matrices that act in flavorspace. These 
matrices are normalized according tr{r a r 5 } = Nf 5 ab and coincide with the Pauli 
matrices in the case of a SU(2) flavor group. The effective quark interaction in the 
soft momentum limit gives rise to the contact interaction ([171) . Let us concentrate 
on the first term of the right-hand side, i.e. 



(D.3) 



Lorentz indices are not explicitly presented. Replacing T A x T A in flDJ) with the 



help of ( p.l|) , the second expression on the right hand-side of (|D . 1| ) gives rise to the 
quark interaction in the vector channel, i.e. 



2N 



qi F i c Yq qi F i c Yq 



Rearranging the product of the 7 matrices, i.e. 



Til 



~ ll k l im + (75)^ (75) 
1 „ „ 1 



(D.4) 



(D.5) 



likl 



im ' 



and using (p.2|) , the first expression on the right hand side of ( |D.1|) gives rise to an 
additional interaction in the vector channel. Taking into account a factor (—1) which 
arises from the anti-commuting Grassmann fields q and q, the vector interaction is 
given by 

' J 1 ' ql F l c rqql F lcYq- (D.6) 



2N f 2 



2N 
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One easily verifies that ( p.6|) also holds for the one flavor case Nf = 1. We observe 
that the interaction in this channel changes sign if the number of flavors is increased 
from one to two: It is repulsive for one flavor but becomes attractive for two or 
more flavors if one takes three colors. It is interesting to note that it is repulsive for 
any number of flavors in the large N limit. Since the large N limit is believed to 
give a qualitatively correct description of hadron interactions, it is in this limit that 
the known repulsion between two nucleons in the isoscalar vector channel (i.e., the 
uj exchange channel) could be understood from the point of view of QCD variables. 
In our work we are exploiting the fact that the gluonic four-quark interaction is 
attractive for Nf = 2 and N = 3. But as stressed this is not the whole story 
because of the possible contribution coming from the collective degrees of freedom 
as discussed in ||14|| . Even if the gluonic four-quark interaction were repulsive to 



start with, as long as it is not too repulsive, it could be overpowered at some density 
by the attraction due to the "dropping" Vq mode as the chiral phase transition is 
approached. This would be the scenario if the chiral transition mechanism in terms 



of hadronic variables discussed in |14j were realized 



The first expression on the right hand side of ( |D.1| ) also induces interactions in the 
scalar and the pionic channel, respectively, i.e. 

gl F l c g gl F l c g - -^-gl c 7 5 r a g gl c 7 5 r a g . (D.7) 



2Nf r r ^ 2N f 

Note that this sign of the scalar interaction is crucial for quark condensation. 

In the following, we specialize to the two flavor case, i.e. Nf = 2. The gluon induced 
interaction that involves charge conjugated quark fields is 

W 6 rl (T*f Y (rf ql 5ik T* Y 4 • (D.8) 

We wish to estimate the strength this contact interaction assigns to the quark in- 
teraction ([D in the (axial) vector diquark channel. For this purpose, first note 
that 



T A A 1 

ap T 7 <5 = T Pa T ~,6 = ^CAaStA-yP + ••• , (D.9) 

Irllik = -jTflTrt - g Wit)* ('fldrk ( D ' 10 ) 

- hllrk + (7 5 )« (7s) rfc + ••• , (D.ll) 



where we used ( p.lp to derive ( |D.9| ). The (axial) vector diquark composite field 



which we investigated in section |] is completely anti-symmetric in flavor space. 
Noting that = —iTfir 2 and using (p.2j) , one finds 

S r iS ik = --eue rk + ... . (D.12) 
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Combining (|D.9| , p.lOj , p,12| ) and taking into account the additional factor (— 1) 



from the interchange of Grassmann fields, we finally obtain the desired result 

- &)? e ik e^ (75) gf q?e lm e^r (Tfe) (D-13) 



1 

4/V 



+ 



D.2 The Fierz rearrangements and approximations 



There are a variety of Fierz rearrangements that can be made on eq. flD]| that could 



lead to different results depending on approximations made in the calculations. The 
expressions given above (eqs. (|D.6| ) and (P.7|) ) represent one such rearrangement 



and are exact as they stand. Note that they involve no color generators and hence 
can in certain approximations be completely saturated by color singlet meson fields. 
An alternative rearrangement was used by Vogel et al |H| wherein even the second 



term of (|D.1| ) is "Fierzed," thereby obtaining an expression that contains terms 
involving the color generators. This is also exact. Thus in the operator form, it is 
identical to our Fierzed form. However the complete saturation of this form in the 
same approximation as with our expressions must require colored states in contrast 
to our formulas. 

An exact field theoretic treatment of the four-Fermi interaction ( p.3| ) should of 



course yield the same result whether one uses our expressions or those of ref.pl 
However if one does approximations in treating this effective four-Fermi interaction, 
there is no reason why one should expect the same result. Indeed if one does the 
conventional mean-field approximation on the two forms, one sees immediately that 
the results would be different, even qualitatively in some cases. Thus a comparison of 
the operators of the same form (e.g., same quantum numbers) would be completely 
meaningless even though in perturbation theory, one should get the same result 
(e.g., the quark propagator) order by order. In the next subsection, we describe the 
strategy - or approximations - that we adopt in handling the effective four-Fermi 
interaction that arises when gluon-mediated interactions flow to the Fermi surface. 
We admit that we have no rigorous proof that our approximation is the best for the 
problem but we believe that it catches the essence of the physics we are developing 
in this paper. 

D.3 The approximation of weakly coupled mesons 

Since the four-Fermi interaction ( p.l| ) is an effective local Lagrangian, we would 
like to work with physical meson fields and approach the phase transition in terms 
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Fock 



Hartree 



+ 



(mean field level) 



Q ab 

a b 



+ 



+ 



(meson loops included) 



Figure 5: Panel a: The leading 
order perturbative diagrams con- 
tributing to the quark propaga- 
tor (solid line) in bosonized quark 
theory; Panel b: leading and 
next-to-leading order of the semi- 
classical expansion of bosonized 
quark theory. 



(a) 



(b) 



of hadronic variables. For this, we find it convenient to set up an approximation 
scheme by introducing auxiliary fields representing quark anti-quark collective fields 
(bosonization). In order to make our argument as general as possible, consider the 
general case of the nonlocal form 



Si ■— q-i ^fkQk G AB qi Tf m q, 



(D.14) 



where k is a master index representing space-time as well as internal degrees of 
freedom. Calculating e.g. the quark self energy in perturbation theory with respect 
to G AB , the s- and t-channel contractions contribute with equal weight and give rise 
(at leading order) to the Hartree and Fock contribution as depicted in figure |5]a. 

Now such nonperturbative phenomena as quark condensation or dynamical mass 
generation take place in a strong-coupling regime which cannot be accessed by per- 
turbation theory. To describe such phenomena, it is more advantageous to rewrite 
the interaction ( [D. 14| ) by means of auxiliary meson fields, i.e. 



exp{-S/} = J VT, exp{-S , A /} 



-i 



Am 



i qiq m £ 



S m — — 

2 V Jik,lm 

The complete partition function can then be written in the form 



1 711 



Z = J VqVqVY, exp j-g, (jS x ) , fc q k - S M } 



(D.15) 

(D.16) 
(D.17) 

(D.18) 
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where Sq 1 = ifi—m is the tree-level inverse quark propagator. Non-perturbative ef- 
fects can be addressed, by integrating out the quark fields and treating the functional 
integral over the "meson" field S in a semi-classical approximation. Expanding the 
action of the effective meson theory up to second order in the meson field S = a + Eo 
yields 

Z = jva exp{-S e// } , (D.19) 
Seff = \<Jik (M- 1 + n) <7 {m - S ml a lm - lnS" 1 + O (a 3 ) , (D.20) 
where S m i is the "constituent" quark propagator with S as self-energy, i.e. 

S- 1 = ip - m + iS , (D.21) 

and 

Rik,lm = Ski S m i ■ (D.22) 

If one chooses in the bosonized action So in such a way that the linear terms in a 
vanish, the quark self-energy will be given by a ladder summation of the Hartree 
terms as shown in the upper graph of figure [5|b. The Fock terms of the perturbative 
expansion are generated, if Gaussian fluctuations o around the mean field S are 
included as in the lower graph of figure [5]b. 

There are several options available depending on how one bosonizes the interaction 



1. One may choose to bosonize eq.( |D.3[ ) directly and then do a semi-classical 



approximation in the mesonic functional integral. This procedure would how- 
ever be a poor approximation for the physics we are interested in. The reason 
is that due to vanishing Hartree terms in this case, at the classical level, it 
would yield a trivial result for the quark self-energy, whereas the nontrivial 
phenomena (such as quark condensation, constituent quark mass generation 
etc.) would occur at semi-classical level or higher. 

2. Suppose one Fierz transforms (partially), thereby rearranging the interaction 
terms (partially). In this case, the role of the Fock and Hartree terms involved 
will be (partially) interchanged. Subsequent bosonization and semi-classical 
approximation of the mesonic functional integral would correspond to a par- 
ticular summation of selected perturbative diagrams. Now the crucial point 
is that different choices of the Fierz rearrangement would differ in the subset 
of Fock and Hartree diagrams which are summed up. Such different choices 
would correspond to doing different physics. The examples are the choice 
adopted in this paper and that adopted in . 
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The first case to consider is the complete Fierz rearrangement of the 
interaction ( p.3| ) studied in |TT|. This choice would reproduce the impor- 



tant Fock terms at tree level. Doing (weak coupling) perturbation theory 
would bring minor corrections from the meson fluctuations. Although 
the semi-classical expansion induced by this choice of Fierz rearrange- 
ment may be rapidly converging in weak-coupling regime, this must not 
be the case for the coupling strength where quark condensation effects 
are expected to come into play. Indeed, an evidence for this shortcom- 



ing can be found in [18]: the Fierz rearrangement of [31] implies a ratio 
of the interaction strengths of scalar and vector channel, respectively, 
of G s /G v = 2 whereas a ratio of G s /G v <C 1 seems to be needed to 



reproduce the ir and p meson masses [|18] 



Consider next the Fierz rearrangement adopted in this paper. We argue 
that doing mean-field approximation with this form corresponds to doing 
the large N c , Nf approximation, i.e. 

N f = uN c , v = fixed const , N c -> oo . (D.23) 

In this limit, the effect of meson fluctuations will be suppressed^. Note 
that the limit ( P-23| ) does not spoil the property of asymptotic freedom 
of Yang- Mills theory for v < 11/2 (note that one has v = 1 for QCD 
with three quark flavors). In the particular case (|D.23| ), the four quark 



interaction strength is solely proportional to 1/N C . Bosonization of these 
interaction terms gives rise to large meson masses, i.e. proportional to 
N c . Mesonic kinetic terms are produced by the logarithm of the quark de- 
terminant. These kinetic terms therefore also acquire a large pre-factor of 
order iV c . These considerations tell us that, in our case of the Fierz rear- 
rangement, meson fluctuations are suppressed by a factor of 1/N C thereby 
justifying our semi-classical expansion in the particular limit ( |D.23| ). 

For practical applications (N c = 2), it is not clear whether a truncation 
of the large N c , Nf expansion (p.23|) at classical level is a truly good ap- 
proximation. In particular at zero density, a higher-order approximation 
may be needed to generate the masses for the vector mesons. Rather than 
justifying the large N c , Nf expansion for the case of interest, e.g. N c = 3, 
Nf = 2, we here adopt a more coarse-grained point of view and argue 
that in view of the results in |[l4j] (see discussion in section |5|) it seems 



worthwhile from a phenomenological point of view to introduce in the 



4 



To those who are familiar with Walecka mean field theory of nuclear matter, this procedure 
is easily comprehensible. One could interpret this in terms of Landau's Fermi-liquid fixed point 
theory with the instability associated with the ISB and color superconductivity interpreted as the 
break-down of the Fermi-liquid ground state. See ]32] for a discussion on this matter. 
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case of finite densities the mean-field Vq and leave the strength of the Vq 
fluctuations as a variable parameter. The goal of our paper is to explore 
the behavior of the finite density quark matter by varying the strength 
in this particular channel rather than relying on the strength provided 
by the effective four-Fermi interaction of one-gluon exchange quantum 
numbers ( p.3| ). 



E Fermion determinants 

Let us first quote two useful results, which will be used below, i.e. 

2 



det D (a# + (3) 
det D (e 7 W + ajT + (3 



a 2 k 2 + (3 2 
(efc+ 



ak ) 2 



+ (e + a)V + (3 



(E.l) 
(E.2) 



where a, f3, e are constants, and := (k^,k) T , k^ = k ± iv. The subscript D 
indicates that the determinant is over Dirac indices. The equations (|E.1| , |E.2|) can 
be readily checked e.g. by resorting to the explicit representation of the 7-matrices 

(ID- 



The determinants of interest possess the structure 



T : 



det 



cDCF 



(E.3) 



where the building blocks of the 2x2 matrix in the q and q c space carry Dirac, 
flavor and color indices. The subscript cDCF indicates the determinant extends 
over either space, i.e. c charge conjugated space, D space of Dirac indices, C color 
and F flavor space. Exploiting the simple structure of the diagonal entries, T can 
be reduced to 
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kl 
ki + a 2 



12 



det 



DCF 



12 



det 



DCF 
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a 



k A (Jk + - a 
K 2 5 2 



ki 
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At 

j/(r + ( x) rt 



(E.4) 



where we have used A = 5T. Let us specify the color, Dirac and flavor structure of 
the diquark vertex, i.e. 

= ^e iklD , (E.5) 

where a, (3 = 1 ... 3 are color indices, and i, k = 1, 2 are flavor indices, respectively. 
Dirac indices are not explicitly shown. The Dirac structure will be specified below. 
Since the propagator (^ + — a) -1 is diagonal in color and flavor space, one obtains 



r # 



a 



-1 



F f 



Id 



(r 



— a 



Y D x l F x P ( 
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(E.6) 
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where lp is the unit matrix in 2 x 2 flavor space, and Pq = diag(l, 1, 0) is a projector 
in color space. We therefore obtain 



T 



2 , a 2 



K 
det^ 



12 



$~ - a + 



k 2 _ + a 2 ' 

k 2 5 2 



k 2 + + a 2 

Investigating the various types diquark composites, i.e 

75 



7d (> + + a) 7d 



Id 
Id 
Id 
Id 



Id 





7°7s 



7 



(scalar) 
(pseudo-scalar) 
(vector) 
(axial- vector), 



(E.7) 



it is straightforward with the help of ( |E.1| , [E.2| ) to calculate the final result for the 
fermion determinants @7|), (glj). 
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